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O We have reexamined the electromagnetic corrections to the weak inter- 

action rates for the transformation of neutrons to protons, and protons to 
neutrons, in the early universe, before freeze-out. We derive compact ex- 
pressions for these rates in terms of thermal expectation values of products of 
fields, and we give explicit constructions of the terms to order e 2 . We disagree 
in several respects with results in the literature. 
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I. INTRODUCTION 



In the early universe, in the era just before nucleosynthesis, the weak interaction rates 
combined with the expansion rate determine the freeze-out neutron-proton ratio at a tem- 
perature on the order of 0.3 MeV. This freeze-out ratio in turn essentially determines the 
primordial He abundance.^] In the region of temperature 0.3 — 3.0 MeV, in which the weak 
interaction rates are not fast enough to preserve chemical equilibrium, it is necessary to cal- 
culate these rates reasonably accurately. The Born approximation rates, as used by Peebles 
H in the first calculations of the freeze-out ratio and in the later refinements by Wagoner 
||, nearly suffice for the purposes of this calculation. But it is possible that the primordial 
He abundance will be measured to such a precision that the electromagnetic corrections to 
the rates will play a role in the interpretation. 

There have been a number of calculations of the electromagnetic corrections to the freeze- 
out ratio Q 0. These corrections come from both Coulomb and radiative effects as 
modified by the hot plasma in the early universe and depend on its temperature. Roughly 
speaking, the method that has been used in the calculations has been the calculation of 
all of the T matrix elements in which n — > p or p — > n, as calculated to order e 2 Gy V . The 
squares of these T matrix elements are then weighted with the appropriate statistical factors 
for the incoming and outgoing electrons, positrons, neutrinos, anti-neutrinos and photons, 
and integrated over the phase space for the particles. We call this the "exclusive channel" 
approach. 

In the present paper we re-address the electromagnetic correction problem. Our motiva- 
tions are several: 

1) There is unresolved disagreement in the literature over the correct treatment of one ar- 
tifact of the exclusive channel approach, the so-called "temperature-dependent wave-function 



renormalization" || [[17]]. And even if the differences with respect to the the electron's 
wave function renormalization were resolved, there is, in exclusive calculations done in Feyn- 
man gauge, the necessity of considering the proton's "temperature-dependent wave-function 
renormalization" as well. This has not been dealt with in the calculations in the literature, 
which were all done in Feynman gauge. 

2) In the exclusive approach followed by all previous authors, two reactions have been 
left out that give contributions of the same order as the corrections that have been included, 
namely the reactions u + e + -\-n^->-p + 'j. We believe that this is mere oversight, originating 
from the fact that these reactions are not corrections to processes that take place in the 
absence of the photon coupling, with a photon added. 

3) It would appear that the reported results of the two most complete reworking of 
the electromagnetic effects, [BJ, J7J, do not obey the appropriate detailed balance relation 
between proton and neutron rates. 

4) The rate correction problem is completely well defined within the standard model 
of the weak and electromagnetic interactions, and it is addressable through perturbation 
theory. Indeed, for the temperature-dependent corrections considered in the present paper, 
the weak interaction can be taken to be of the local four-Fermi form. We shall endeavor to 



1 For a simple pedagogical treatment of the abundance calculation see Ref. (T| 
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present a definitive and final result. In view of all of the past confusion, we shall take pains 
to give complete arguments, to give the analytical results in detail, and to explain all of the 
steps. 

We shall develop the perturbation theory for the neutron and proton disappearance rates 
directly, without going through the intermediary of the calculation of T matrix elements for 
exclusive reactions. Instead, we begin from equations for the inclusive neutron and proton 
transformation rates, expressed as the appropriate integrals over thermal expectation values 
of electron field operators. We then apply the standard methods of thermal field theory 
to obtain the electromagnetic corrections of order e 2 and to reduce the answers to a form 
allowing easy computation. We prove that the simple detailed balance relation between the 
proton and neutron rates remains valid to order e 2 . This relation serves as a useful check on 
results. There is no role for "temperature-dependent wave-function renormalization" in our 
approach. Our calculation ends with numerical results that are quite different from those of 
previous treatments. 

We shall assume complete thermal equilibrium for each species. The effects of a non- 
thermal neutrino distribution have been dealt with separately elsewhere || , fTfl . Their effect 
on the terms of order e 2 , a correction to a correction, is not consequential, particularly so 
since the deviations from a thermal spectrum become appreciable only at temperatures very 
near the freeze-out point. 

However the system can be arbitrarily far from chemical equilibrium. It will suffice to 
consider a single neutron in the medium in order to calculate the neutron transformation 
rate, and a single proton in the medium to calculate the proton transformation rate. The 
electron and neutrino chemical potentials are taken to vanish. 

II. GENERAL FORMULATION 

The rates of neutron and proton appearance and disappearance can be expressed in terms 
of the inclusive rates of v e and v e appearance and disappearance. Since we are concerned 
only with the modifications of the rates which arise from plasma effects that involve low 
energies — which are 'soft processes' — we may take the nucleon mass to be infinite] 2 ]. We 
take the proton or neutron to be situated at the origin. We begin by taking the initial 
nucleon to be a neutron. As described in detail in Appendix A, the rate for a neutron to 
change into a proton may be expressed as 

Vn= Wf + ^) loo ^ El n(Eu) Loo ^ e+tEVt W(t) ' (2 ' 1} 

where 

w(t) = (Mo,t)T4t)4(o)T + (o) ) t . (2.2) 



2 The zero temperature rates, which we shall delete from our results, involve a domain of very 
hard virtual photons that requires the inclusion of nucleon recoil. The zero temperature radiative 
corrections are summarized in Ref's. || and jj]]. There are also small temperature-dependent recoil 
effects that do not involve the electromagnetic couplings. These are evaluated in ref. ||. 
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Here the angular brackets with the T subscript denote the thermal average including all 
plasma and electromagnetic radiative interactions as well as the single particle state of 
the (infinitely heavy) neutron placed at r = 0. The electron field operator is denoted by 
ip e {x) and T±(t) are the isospin raising and lowering operators, with all the operators in 
the Heisenberg picture — their time dependence is controlled by the full Hamiltonian of the 
system which also defines the thermal average. We shall assume that the chemical potentials 
of all the leptons vanish so that they involve the Fermi occupation density^ 

Note that 

n(-E) = e pE n(E) = 1 - n(E) (2.4) 

is the Pauli blocking factor for a particle produced with (positive) energy E. 

The integral over positive neutrino energies E v in the rate T n describes neutrinos ab- 
sorbed from the thermal bath with a population governed by n(E u ). Including corrections 
of order e 2 , it is the sum of the rates for the processes: 

1. v + n — > p + e~ 

2. u + n + j — > p + e~ 

3. z/ + n^p + e~ + 7 

4. v + n + e + — >p + 7 

The last process in this list, u + n + e + — > p + 7, and the reversed process, have been omitted 
in all previous calculations in the early universe application. Their contribution will turn 
out to be at about the same level as some of the processes considered in Ref's. @] - @. 

The integral over negative neutrino energies E v in the rate T n describes antineutrinos 
emitted into the thermal bath with the Pauli blocking factor n(—\E v \). Again to order e 2 , 
it is the sum of the rates for the processes: 

1. n —>■ p + e~ + v 

2. n + 7— > p + e~ + v 

3. ?2^p + e _ + 7 + z/ 

4. n + e + — ► p + v 

5. n + e + + 7^p + z/ 



3 These apply to the free neutrinos and to the electrons once the electromagnetic corrections 
appear as a perturbation. The electron distribution is, of course, modified by the electromagnetic 
interactions, but this is accounted for automatically in our formalism. 
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6. n + e + ^p + / y + u 

Again as shown in Appendix A, the rate for the disappearance of an initial proton is 
given by 

T p = 0j~2 (9v + dE v El n(E u ) dt V(t) , (2.5) 

where 

V(t) = (^mT + (0)MO,t)T4t)) T . (2.6) 

The integral over positive E u describes the rates involving the absorption of antineutrino, 
the processes in the second list above with the direction of the reaction arrow reversed. The 
integral over negative E v describes the emission of neutrinos, the reverse of the first list of 
processes above. For future reference, we note that the simple change t — > — t gives 

Tp = (# + ^) loo ^ El H{Eu) J-oo ^ ^ Vit) ' (2 - 7) 

where V(t) = V(—t). We make use of the time-translation invariance of the expectation 
value to add the time t to all operators to obtain the form 

V(t) = (4(0, t)T + (t)M0)T-(0)) T ■ (2.8) 

This expresses the proton rate T p in exactly the same form as the neutron rate T n except 
for the interchanges ipl <-> tp e and T + <-> T_. 

It should be emphasized that our formalism takes account of all these individual pro- 
cesses and unifies them in terms of a simple expression. This expression forms the basis in 
which the plasma effects on the radiative corrections are easily computed in a complete and 
unambiguous fashion. 

To the order e 2 to which we work, the neutron and proton proton rates, T n and T p , obey 
the detailed balance relation 

r p = e^ A r n , (2.9) 

where 

A = M n — M p (2.10) 

is the neutron-proton mass difference. Thus only one of the two rates T n or T p need be 
calculated. 

To prove the detailed balance statement, we first define the projection operators P n and 
P p for the point-like heavy neutron and proton states located at r = 0. In terms of these 
operators, we have the explicit expressions 

W{t) = Z- 1 Tre-^e(O,i)T_(i)^(0)T + (0)P n , (2.11) 

and 
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V(t) = Z p l Tre-P H 4(0)T + (0)MO,t)T4t)P P , (2.12) 

where 

Z n = Tr e-? H P n , Z p = Tr e^ H P p . (2.13) 

The projection operators P n and P p are time independent (commute with H), and they 
commute with the electron fields, while 

T_(t)P p = P n T_(t). (2.14) 

Using this result, the cyclic symmetry of the trace, and 

ip e (0,t)T_(t)e-P H = e- pH i) e {0,t - i(3)T_(t - i/3) , (2.15) 

we obtain 

V(t) = Z p 1 Tr e"^i(0, f - - ^^(0)^(0) P„ 

= (f?) W(t-i(3). (2.16) 



Placing this result in the rate formula ( |2.5| ) for the proton and shifting the time integration 
variable t — > t + i/3, we find that 



r _ G w 



v Z n (2tt) 2 



9v + 39a) / dE v E 2 v n{E v )e pE " / tfte - *^*^). (2.17) 



Finally we recall Eq. ( |2.4| ) which states that n(E u ) exp{(3 E u } = n(—E u ). Accordingly, we 
reflect the neutrino energy integration variable, E u — > — to place this result in precisely 
the form ( |2.1| ) of the neutron rate except for the initial overall factor, and hence 

r P =(^Vn. (2-18) 



In lowest order, 



^ = e^ A , (2.19) 

Zj p 



and the detailed balance statement Q2.9| ) is obtained to order e 2 if there are no corrections 
to this ratio to order e 2 . It is a familiar result from statistical mechanics that there is no 
order e 2 correction to the partition function for a dilute gas. The same considerations show 
that there is no e 2 correction when a single charged particle is introduced into a plasma. 
This is made explicit in Ref. [2^] where it is also shown that the first correction is given by 

bZ n 1 . 
-^Z = -p aKD , (2.20) 

where kd is the Debye wave number. The nature of this result is clear: The introduction of 
a charged particle into a plasma alters its electrostatic field from a pure Coulomb potential 
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e/47rr to the Debye screened potential (e/Anr) exp{— n^r}, and thus changes the particle's 
self-energy by 



5E = lim - - 

r^o 2 r 



-K D r 



- 1 



a Kn 

2 



(2.21) 



giving SZ p /Z p = —(3 5E. For the electron-positron plasma that is of our concern, the squared 
Debye wave number is given by 



K 2 D = (3e 2 Af = e 2 



dp, 



(2.22) 



iu=0 



where I is the leptonic charge density, 



(2tt) 3 



1 



exp{(3[E(p) -//]} + ! exp{/3[E(p) +//]} + ! 



^2.23) 



In the high temperature limit in which the temperature is much greater than the electron 
mass, T>m, this gives k 2 d = 4naT 2 /3 and so the plasma correction (12.20Q is a negligible 
change 



5Z„ 



a 



3/2 



0.00064 



^2.24) 



At lower temperatures, this correction vanishes exponentially (~ exp{— /3m}) as the tem- 
perature becomes small in comparison with the electron mass, and is thus even smaller. 

To demonstrate the workings of our formalism, we exhibit the rates in the absence of 
plasma interactions and radiative corrections. In this limit, 



T_(f) 



-iAt 



T_(0) 



^2.25) 



and using the free field correlators ( |A17| ) for the thermal values of the electron fields in 
Eq. (^T]) we get 



ri 0) = G 2 W (g 2 v + 3g 2 A ) \ J ^ [n(-E) (E - A) 2 n(E - A) 



(2 

+n(E) (E + A) 2 n(—E — A) 



(2.26) 



where we use the short-hand notation E = E(p). The first term in the square brackets in 
the region m e < E < A describes free neutron decay n —>■ p + e" + v with the final leptons 
having the Pauli blocking factors n(—E), n{— |A — E\). For E > m e , the term corresponds 
to v + n — > p + e~ with the initial thermal neutrino population governed by n(E — A) 
and the final electron having the Pauli blocking factor n(—E). The second term in the 
square brackets represents the process e + + n — > p + v with the initial positron drawn from 
the thermal bath with population factor n(E) and the final antineutrino having the Pauli 
blocking factor n(—E — A). In this limit of no corrections, the proton disappearance rate 
( p.5[ ) becomes 
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r£ 0) = G 2 W (gl + 3g 2 A ) I / |g| [n(£) (£ - A) 2 n(A - E) 

+n(-E) (E + A) 2 n(E + A)] , (2.27) 

with the terms in the square bracket describing the previous neutron processes in the reverse 
direction. Since 

n(±E) = e T ? E n(TE) , n(A ^ E) — e ±E E~ /3A n(±E — A) , (2.28) 

these rates do indeed obey the detailed balance relation Fj® = exp{(3A}F^ . Note that the 
neutron and proton rates F n and F p are related by the interchange^ A <-> — A. 

In view of the detailed balance relation (|2.9|) , we shall concentrate on the neutron disap- 
pearance rate F n . To write the lowest-order result in a compact form, we define 

X (E) = n(E - A)n(-E)(E - A) 2 , (2.29) 

and integrate over the electron solid angle to obtain 

^ = 4 (# + /o ^ ^ + • (2 ' 30) 



The general from for the neutron disappearance rate is given by Eq's. Q and (|J 



Since the rate involves the operators ?/> e (0, t)T_(i) and ^](0)T+(0) which are electrically 
neutral (they commute with the electromagnetic current), the rate is gauge invariant. Since 
we treat the nucleons as being infinitely heavy, it is convenient to perform calculations in 
the radiation or Coulomb gauge, and this we shall do. In this gauge and with a very heavy 
proton, the second-order electromagnetic correction involving the proton alone is simply 
a mass shift that is removed by the usual renormalization. To the order e 2 to which we 
work, and with the radiation gauge, the only non-trivial corrections come from the Coulomb 
interaction between the electron and proton and the radiative and Coulomb corrections to 
the electron (and positron) system. 

We define the T = rates as the rates that one would get summing the contributions 
from all of the processes listed in the introduction, where thermal distributions are taken for 
all initial and final electrons, positrons, neutrinos and antineutrinos, but the temperature 
is taken the vanish everywhere else. This completely excludes participation by the thermal 
bath of photons. Thus the only real photon effect in our definition of the T = rates is 
a process with a bremsstrahlung photon in a final state. The T = terms include all the 
radiative and Coulomb corrections that produce ultraviolet divergences which are removed 
by the usual renormalization of the vacuum amplitudes. 



4 This equivalence for the unperturbed rates follows directly from the alternative form (^]) for 
T p . Because of charge-conjugation invariance, the free-field expectation values of ip^ (t)ip(0) and 
if}(t)^ (0) are equal (in the Majorana representation of the Dirac matrices). Replacing T-(t) by 
T + (t) changes the factor e +tAt to e~ tAt . The expectation value of T_T + in the neutron state equals 
that of T + T_ in the proton state. 
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We shall subtract all such "T = 0" corrections to define and work only with temperature- 
dependent corrections to the rates that are free of the renormalizations needed to deal with 
the ultraviolet infinities. 

This definition of T = corrections agrees with the definition as used in Ref's. Q 
|7j except in one regard: If we consider that part of the reaction rate for the reaction 
n + e + + i>—>-p + '-fm which the Bose factor for emission [1 — exp(— /3u;)] _1 is replaced by 
unity, this term falls under the above definition of a T=0 term, but it is not included in the 
T = terms of Ref's. |4[] 0, for the reason that the process was not considered at all in 
these works. 



III. ELECTRON-PROTON INTERACTION 

We shall first look at the plasma correction to the Coulomb interaction between the 
electron and proton, 

H c (t) = -J (dV)(d 3 r) *l(v', t)%(v', t) f_ V&r, flVfefo t) . (3.1) 

We start with this piece because it serves as an easy introduction to the method that will 
later be applied in the more complex electron self-energy correction. To account for the 
Coulomb perturbation in first order, we temporarily continue to imaginary time, t — > —it, 
with t > 0. Then the perturbative rules of thermal field theory in imaginary time apply, 
and the electron-proton Coulomb interaction correction to Eq. (|2.2j) becomes 



(-it) = - | o /3 rfr^(^(0,-^)T_(-^) J ff c (-zr')^ t (0)r + (0)) + ^ , (3.2) 



where the (•■•)+ denotes time ordering in the imaginary time. The operator Hc{—iT') makes 
no contribution when r' > r because the ordering places it adjacent to the neutron state 
on which it vanishes. Taking this into account, we may now continue back to real time and 
write 

W ep (t) = -i J*dt'(MO,t)T4t)H c {t')4{0)T + (0)) T , (3.3) 

where now no time ordering is involved. Since the perturbation is explicitly taken into 
account, we may now use the non- interacting results. Since the expectation value of 
^p(r', t)ty p (r', t) is sharply localized at r' = 0, the coordinate r' may be taken to vanish 
in the Coulomb potential factor. Hence the total charge operator of the proton appears 
which acts on a proton state of unit charge, and within the expectation value we have 

T_(t) / (d 3 v')^ p (v',t)%(v',t) T+(0) =T_(t)T + (0) - e +iAi . (3.4) 

Moreover, the electron field thermal expectation value may be replaced by the product of 
free fields defined in Eq. (|A13| ) of Appendix A. Thus 



W ep (t) = +ie +iAt C dt' /'(d 3 r)-^- r trS , W(-r,t-t / )7 5 , W(r,t , )7 : ( 3 - 5 ) 
Jo J Aix\y\ 
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and using the explicit form (|A17|) for the functions that appear here, 



W ep {t) = +4e 



(d 3 p) 1 f(d 3 p') 1 



(2tt) 3 2E J (2tt) 3 2E' (p - p') 2 
EE' + m 2 + p • p' 



E-E> 

-n(E)n(E' 

EE' - m 2 - p • p' 

" E + E' 



n(—E)n(—E') 



e - iE 't _ e -iEt 



e iE't _ e iEt 



n(E)n(-E') 



e iEt _ e -iE't 



(3.6) 



where E = E(p) and E' = E(p'). 

In placing this result in formula ( |2.1| ) for T n , the time integration produce energy- 
conserving 5 functions. When the resulting neutrino occupancy factors appear in the form 
of n(E — A) multiplying n(—E) or n(—E — A) multiplying n(E) [or the same forms with 
E — > E'], we have a "T = 0" structure that, as described above, is to removed. These are 
just the neutrino and electron occupancy factors that appear in the free rate ( |2.30| ). The 
terms involve a further factor of n(—E) = 1 — n(E) [or n(-E') = 1 — n(E')} with the "1" 
parts giving the (Coulomb corrected) "T = 0" part that is to be removed. Our definition 
thus produces the temperature-dependent rate correction 

2 r(d 3 p) 1 f(d 3 p') 1 e 2 



ra, = G w (9v + 3$ - J 



(2tt) 3 2E J (2tt) 3 2E' (p - p') 



EE' + m 2 + p ■ p' 



n{E)n(-E') (E' - A) 2 n(E' - A) 



E — E' 

-n{-E)n{E') (E - A) 2 n(E - A) 
+ n{E)n{E') \(E' + A) 2 n{-E' - A) - (E + A) 2 n(-E - A) 
EE' - m 2 - p ■ p' 



-2 



n{E)n(E') (E + A) 2 n(-E - A) 



E + E' 

+n(-E)n(E') (E - A) 2 n(E - A) 



(3.7) 



We make use of the definition (|2.29| ) of x(E)i use pdp = EdE, and perform the integrals 
over angles to obtain 



n ,ep 



e 2 G 2 



W ( 2 



(2tt) 5 



g 2 v + 3s 2 / dE / dE' 



(E + E') 2 In fp + P " 



4pp' 



E'-E 



p — p 

n(E') X (E) - n(E) X (E') + n(E') X (-E) - n(E) x (-E') 
App' — (E' — E) 2 In ' '" ] 



p — p' 



E' + E 
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n(E') X {-E) + X (E) 



(3.8) 



Each term in the integrand has a Fermi distribution function which falls exponentially for 
large E or E', so that there is no ultraviolet divergence. The terms in the second square 
bracket are odd under the interchange E <-> E' and thus vanish when E — E' so as to remove 
the singularity in the overall (E — E')^ 1 factor, leaving an integrand that is integrable at 
E = E'. Thus we are left with a well-defined double integral that may be done numerically. 
However, to save writing, and to place this result in a form that will prove useful later, we 
exploit the symmetry of the double integral to express the result in an asymmetrical form 
with the understanding that the potentially singular terms are defined by a principal part 
prescription: 



n ,ep 



(27T) 



dE 

m Jm 



E 



E' 2 - E 2 



dE'n(E') [ X (E) + X (-E) 
4p P > - (3E< 2 + E 2 ) In j 



(3.9) 



It is not difficult to confirm that the principal part prescription is equivalent to the original 
symmeterization of the double integral. 

Making use of n(E) = e~ l3E n(— E), we can write the first form 03-71) of the result as 



r( T ) 

n ,ep 



G 2 W (g 2 v + 3<£) - 



(d'p) 1 /"(dV) 1 



(2tt) 3 2E 



(2tt) 3 IE' (p - p 



A2 



EE' + m 2 + p • p' 



n{E) n(E') 



E-E' 

{{E' — A) 2 n{—E' + A) - (E — A) 2 n(—E + A)| 
+ \(E' + A) 2 n(-E' A ^ 1 A ^ 2 



EE' 



m 



PP 



E + E' 

+(E — A) 2 n(E — A) | 

In this guise, it follows immediately that the formal change A 



A) - [E + A) 2 n(—E — A) j 
£) n(E'){e /3A (E + A) 2 n(E + A) 



-A produces 



(3.10) 



(3.11) 



This is just the detailed balance relation ( |2.9| ) that relates the neutron and proton rates. 
Since the subtracted "T = 0" terms are proportional to the uncorrected rates that obey this 
relation, we conclude that the corresponding "T ^ 0" part of the Coulomb electron-proton 
correction to the proton rate may be obtained simply by the interchange A — > —A. That 
is: 



r (T) = p(T) / A ^ —A) 
p ,ep n ,ep\ J 



(3.12) 
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IV. ELECTRON SELF-ENERGY AND RADIATIVE CORRECTIONS 



We turn now to examine the remaining leading-order electromagnetic correction that 
involves only the electron (positron) in the neutron rate defined in Eq's. ( [2.1]) , ( |2.2| ). Since 
this effect involves no electromagnetic interaction with the proton, we have, effectively, 

T_(t)T+(0) -> e lAt . (4.1) 

It is convenient to translate the energy integration variable by E v — > E = E v + A to 
re-express the neutron rate formulas ( |2.1| ), fl2.2|) as 

Gl 



r ne _ e = — ^ (<# + 3g\) j_2 dE (E - A) 2 n(E - A)W(£) , (4.2) 
in which 

W(E)= dte iEt (i>(O,t)^(0)) T . (4.3) 



The angular brackets now denote the thermal average only for electrons and photons and, 
to reduce notational clutter, we have removed the e subscript on the electron field operators. 
Note that E > describes outgoing electrons. 

As we have just seen, the effects of interactions are most conveniently dealt with in the 
imaginary time formulation. However, the continuation to real time is now no longer trivial. 
Hence we pause to review very briefly the relationship between the unordered real-time func- 
tion ( |4.3j ) and its counterpart, the time-ordered in imaginary time electron thermal Green's 
function. Since the thermal average (• • -)t is the normalized trace Z~ l Tr exp{— /3H} ■ • •, 
the cyclic symmetry of the trace may be invoked to interchange the order of operators by 
passing then through the evolution operator in imaginary time exp{— f3H}. Therefore 

(^(o, t)^(o) ) T = <y (o) m t + %0)) T . (4.4) 

In term of Fourier transforms using the sign convention e~ lEt , this relates the transforms of 
the two operator ordering by a factor of e l3E . Therefore, 

, , r+oo A TP 

(m t)^(0) ) = / — e~ im n{-E) A(E) , (4.5) 

and 

<y (0 )^(0,t)) T = —e-* Et n(+E)A(E), (4.6) 

where A(E) is the Fourier transform of the thermal average of the anticommutator and 
n(E) is the Fermi distribution function with vanishing chemical potential that we have 
already made much use of. The thermal Green's function, time-ordered in imaginary time, 
are computed by the usual quantum field theory rules, but in Euclidean space, t — > — ir, 
and with Bosonic functions periodic and Fermionic functions anti-periodic in the imaginary 
time interval , f3. Thus the electron's thermal Green's function has the Euclidean 4-th 
momentum component = (2m + l)vrT, and using the representation (|4.5| ), we have 
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(d 3 p) 
(2tt) 3 



G(p,ip±)j = j dTe^ T (tfj(O,t = -zr)^(0)) rj 

» dE' A(E') 
3 2n E' — ip4 



(4.7) 



Continuing back to Minkowski space-time, with p 4 — > — zi£ + e = — + ie), e — > + , and 
taking the imaginary part produces 

A(E) = 2 / £^ Im G(p, £ + ie) 7 ° , (4.8) 



(2tt) 3 

and so, in view of Eq's. (|4.3|) and ( |4.5|) , 

W(E) = 2n(-E) |j^ImtrG(p,£ + *e) 7 °. (4.9) 

The thermal electron Green's function has the structure 

G(p, m) = 7mPm + rn + E(p) , (4.10) 

where 74 = — i 7 °, and we shall compute the self-energy function E(p) to first order in e 2 . 
We separate out the divergent terms that are removed by the renormalization process by 
writing 

E(p) = E(°)(p) + E( T )(p), (4.11) 

in which the first term on the right-hand side is the vacuum function and the second 
term vanishes when the temperature vanishes. The vacuum part gives the renormalization 
terms. We avoid dealing with these renormalization effects by computing only the "intrinsic" 
temperature-dependent corrections, the temperature-dependent correction for the internal 
lines of order e 2 graphs that have no divergences. We continue back to Minkowski space-time 
and write 

1 m — 7 p 



7 p + m E(p) 2 — (E + ie) 2 ' 
where now 7 p = ^kVk — 1°E, and E(p) = VP 2 + rn 1 . To our order 

1 

[E(pf~ E 2 ] 



(4.12) 



Im G(p, E + ie) = — {rn — 7 p) \ — Im E(p, E + ie) 



+ E^ (p, E) Im L ■ - ^ 1 (m " 1P) ■ (4 13) 

[E{p) 2 -(E + ie) 2 } J 

In the second term on the right here, only the temperature-dependent part of the self-energy 
function enters — the vacuum contribution is removed by renormalization. To deal with 
this second term, we note that 

1 Id 1 

Im tt = — — — — Im 



[E{p) 2 - (E + ie) 2 } 2 2EdE E(p) 2 - (E + ie) 2 

— S ^rr [S (E - E(p)) -S(E + E(p))] . (4.14) 

2EdE 2E(p) [ v W!! v yy!n v ; 
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To summarize, we recall the definition Q2.29Q of x(E), and use Eq's. (|4.14j) , ( |4.13| ), and 
^9|) in the rate formula (14.21) to obtain 



r 

1 ne—e 



r 2 



(2 



7T 



-(m - 7P) rrn/ rn ? i2 2 ImS (p ; £ + i e )( m ~ IP) 
[E(py - e 2 ] 

+ [5(E-E(p))-8(E + E(p))] 

7T 8 



2E(p) dE 



E~ l x(E) (m — ^p)^ T \p, E){m — jp) \. (4.15) 



In the last line, the derivative of the delta function has been integrated by parts. When 
this derivative acts upon x{E), the remaining factor is gauge invariant, and the result 
corresponds to a temperature-dependent shift of the electron's energy in the plasma. The 
terms that result when the derivative acts on the remaining factors could be associated with 
some sort of "temperature-dependent wave-function renormalization" , but these terms are 
not gauge invariant, and so they cannot be given any independent physical interpretation. 
The first line, involving the imaginary part of the complete self-energy function describes 
bremsstrahlung processes in the plasma and, since an imaginary part is taken, it has no 
ultraviolet divergences. The first line does, however, contain "T = 0" contributions that 
must be removed. This is discussed in the next section and in more detail in Appendix B. 



It should be emphasized that the electron-electron contribution (|4.15| ) given here, to- 



gether with the electron-proton contribution (|3.9|) given at the end of the previous section, 



give a complete and unified description of all the thermal, plasma effects on the leading 
electromagnetic corrections to the neutron disappearance rate. Some distinct and separate 
physical processes may be identified in the total rate contribution such as the energy shift 
mentioned above which will soon be described in more detail. However, other terms, such 
as those that may be associated with a 'wave function renormalization' contribution, are 
not gauge invariant and have no physical meaning whatsoever. Indeed, the electron-proton 
Coulomb correction computed in the previous section is not gauge invariant and has no 
physical significance by itself. We turn in the next section to assemble our results in terms 
of physically significant parts. We do this so as to obtain formulae that are easily read and 
not so long as to tire the eye, but still have well-defined physical significance. The details 
of the calculation of these results appears in Appendix B. 

Before going into the details of our results in the next section, we note that the alterna- 
tive form (|2.7|) shows that the proton rate may be obtained from the neutron rate by the 
interchanges ip <-> -0 and T + <-> T_. Since r e _ e involves no interaction with the proton, 
charge conjugation invariance holds for the electron field thermal expectation values, and 
so they are not altered by the interchange ip <->• ijy . The interchange T + <-> T_ just changes 
the associated time dependence from e tAt to e~ tAt . Hence the proton rate is related to the 
neutron rate by 

r pe _ e = r ne _ c (A — > -A) . (4.16) 



14 



V. RESULTS 



As we have just noted, the first line of Eq. ( |4.15| ) describes real photon processes. One 
piece involves a factor of the thermal photon distribution function 



/(*) 



1 



(5.1) 



and the remaining part does not involve this function. This remaining part has pieces 
corresponding to the real photon emission processes z/ + n^p + e~+7, u + e + + n^p + r y, 
n — >p + z/ + e~+7, and e + + n — > p + v + 7. As described in Appendix B, these pieces 
correspond to parts of the "T = 0" rate that are easily identified and removed to produce the 
desired T^O rate contribution. The remaining part of the real photon processes described 
by the first line of Eq. Q4.15| ) is the photon thermal bath contribution ( |B30| ) computed in 
Appendix B. This result reads 



p(7,7) _ o 
1 ne— e ^ 



SjK^)/;f /w /> 

[n(-£)x(£ -k) + n(£)x{-£ + k)\ F_ + [n{-£)x{£ + k) + n(£) X (-£ - k)\ F. 



Here 



X(E) = n(E - A) (E - A) 2 
according to Eq. (|B29|) . The functions 

F± = A± kB 

appear in Eq. (|B20| ), with 



A 




and 



2£ 2 + k< 



B = 2S\n 0±A -Aq. 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



where 



q = y/£ 2 — m? 



(5.7) 



Although this set of terms in Eq. 



has no ultraviolet divergences, since f(k) k 1 



k as k — > 0, it is infrared divergent. This contribution is rendered finite in the infrared 
by a part of the contribution contained in the last two lines of Eq. ( |4.15| ) []. It is one of the 



5 The fact that the infrared divergences must cancel, order by order, for these rate calculations 



was proved in [18] 
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parts that arises from all the terms that remain when the energy derivative d/dE does not 
act on x{E)- This piece, which involves an overall factor of the photon distribution function 
f(k), is calculated in Eq. ( |B53| ) of Appendix B to be 



r 7 ' z = 9 

1 ne-e ^ / \ 5 



e 2 G 2 



IV 



(g 2 v + 3g 2 A ) j o -f{k) J m d£A{-2n(-£)m-M£)x(-£)} ■ (5-8) 



(2kY 



Adding Eq. ( |5.8[ ) to the real photon rate ( |5.2| ), and writing our rationalized charge in terms 
of the fine structure constant, a = e 2 /47r, gives 



,(7,/) 



2a G 2 W 



71 (2?r) 



<£ + 3<&) jfy/(*) jTde. 



.4 



-jfe-B 



n(-£){x(£-A;) + x(£ + A;)-2x(£)} 



+n{S) {x(S + k) + x{S - k) - 2 x {-£)} 



n(-£) {x{£ - k) - x{£ + k)} + n{£) {x{-£ + k) - x{-£ - k)} 



(5.9) 



The weight k" 1 f(k) diverges as k~ 2 when k vanishes, k — » 0. The set of terms in the first 
square brackets is of order k 2 for small fc. The terms in the final square brackets are of order 
k for small k, but they are multiplied by an explicit factor of k. Hence, this contribution to 
the neutron rate has no infrared divergence. The k and £ integrals are damped for large k 
and £ by the Bose and Fermi phase-space density functions, and so it is also well behaved 
in the ultraviolet as well as the infrared. Since this is the only part of the rate that is 
proportional to the photon density f(k), it is a well-defined, gauge invariant part of the 
neutron disappearance rate. 

The final pieces that arise when the energy derivative d/dE does not act on x{E) in 
the last two lines in Eq. ( [4.15| ) is presented in Eq. ( |B54| ) in Appendix B. Adding this to the 
modification (|3.9| ) of the neutron rate brought about by the Coulomb interaction between 



the electron and proton and again using a = e 2 /4tt gives 



2 a Gl 



n , ep+ee 



W 



it (2tt) 



(g 2 v + 3g 2 A ) dE\x(E) + X (-E)} dE' n(E' 



E 



E 



/2 



8pp' — 



E' 2 + E 2 



In 



p-p' 



E' 2 + E 2 



In 



[EE' + pp'\ 



[EE' — pp'\ — m 




(5.10) 



Since the only remaining term involves the gauge-invariant energy shift shown in Eq. (|5.11[) 
below, and all the other previous contributions are gauge invariant, we conclude that 
Eq. ( p.l(j| ) is a well-defined, gauge-invariant correction to the neutron disappearance rate. 
Just as we have done in the writing of Eq. ( |3.9| ), it is left implicit that the potentially singular 
terms are defined by a principal part prescription. 
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Finally, there is the term when the energy derivative in the last line of Eq. ( |4.15|) acts 
on x(E). In this case, the electron self-energy function is evaluated on mass shell and 
thus describes a gauge-invariant, temperature-dependent energy shift correction which, after 
performing the angular integrations, has the form 



^2 = 4^(^ + 3^) dE J Q dp 



(27T) 



a 



[5 (E - E(p)) + 5(E + E(p))] AE^(p) — X {E) . 



(5.11) 



As explained in more detail in Appendix B, the energy shift AE^{p) is a shift in the 
position of the pole in the interacting thermal electron propagator — it gives the energy 
shift to create an electron (positron) of momentum p in the plasma relative to the energy 
needed to create the particle in the vacuum. The form of Eq. Q5.11J) is precisely the change 
in the free rate ( |2.30p when the energy-momentum relation of the electron (positron) is 
altered. The energy shifts is computed in Appendix B. The results displayed in Eq's. ( B45|) 
and ( ]BT7| ) give 



i - 



m 
4 qp 



In 



q + p 



(5.12) 



The energy derivative in the energy- shift term (pUll) may be removed by an integration by 
parts using d/dE = (E/p) 8/ dp. Using the explicit form ( |5.12| ) of the energy shift, we find 
that 



— 7^3 (?v + r - dE + X(-E) 

7T (Ztt) 6 v ' Jm p 1 

{ / dkkf(k) + / dE'p'n(E' 

JO Jm 



roo £ 



E' 2 - E 2 



(5.13) 



Again it has been left implicit that integrand in the second double integral in E and E' is 
defined by the principal part prescription. 

The non-trivial plasma effects on the one-loop electromagnetic correction to the total 
neutron disappearance rate is the sum of the results that we have enumerated: 



Arl T ) 



p(7,/) +r 



, r(AB) 
n, ep+ee ' 1 ne-e > 



(5.14) 



in which the successive contributions are given in Eq's. ( |5.9|) , ( |5.10|) , and ( 5.13|) . We should 
emphasize one last time that this result follows from a straight forward evaluation of the 
corrections to the basic rate formulas (|2.1|) and (|2.2|) . We have broken up the complete 
result into smaller pieces in large part so as to have reasonably short expressions that do 
have gauge- invariant, physical meanings. But we have not patched together a set of various 
physical processes. Rather, we have started from a well-defined general formulation and just 
worked out the consequences. 

The second set of terms on the right-hand side of Eq. (|5.14|) involve singular integrands 
that require principal part definitions. To write the result in a less singular form, we perform 
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some algebraic rearrangement^ and partial integrations to express the sum of these two terms 
in the form 



r r (AE) _ 2 a G 2 W 

1 n,ep+ee 1 ne-e - ^ ^ xg 



9v + 3g\ 



dE 



X(E)+ X (-E) 



with 



F(E,E') = --In 2 



' p + p r 



p — p 

2 



n'{E>) P -^[E + E'} + n(E') ^ 
p E' pp' 



(5.15) 



E' + 



m 2 E 
E' 2 



+ ln 



( p + p' 
\p-p> 

+ n(E') 



n'(E' 



Em 



p 2 E 



^ (E' 2 +2p 2 +m' 



P 

E 2 + E' 2 _ E 2 ^ 
E + E' pp 1 K ' J 



- n(E') \ AE - + 2E' L(E, E') 
I p 



(5.16) 



Here in the first equality we have explicitly evaluated the simple integration over the photon 
distribution f{k) in Eq. ( |5.13| ), and in the second equality we have used the definitions 



n\E') 



dn(E') 
dE 1 



and 



L(E, E') = In 



/ EE' +pp' + m 2 " 
\ EE' — pp' + m? 



(5.17) 



(5.18) 



This non-singular expression of these rates is a convenient form to use in numerical evalua- 
tions as we shall do below. 

It was shown in Eq. ( |3.12|) that the changes in the neutron and proton rates brought about 
by the electron-proton Coulomb interaction are related by the formal substitution A <-> —A. 
The same substitution relates the total electron self-energy and radiative corrections to rates 
as shown in Eq. (|4.16|) . A glance at the details of the "T = 0" subtractions made to the 
electron self-energy and radiative corrections shows that those for the proton and neutron 
are also related by the interchange A <-> — A. We therefore conclude that the complete 
"T ^ 0" corrections to the proton rate are given by 



Arf) = Arf)(A^-A). 



(5.19) 



'An ingredient needed for this is the identity 

[EE' + pp'} 2 - m 4 _ / EE' +pp' + m 2 \ 2 fp + p'\ 2 
[EE' - pp'} 2 — m 4 ~ \ EE' — pp' + m 2 ) \p - p' J 
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As was discussed in the Introduction, the process z/ + e + + n— > p + 7 has been hitherto 
inadvertently omitted in the literature. This contribution to the neutron rate is identified 
in Eq. ( |B31| ), which we repeat here: 

T nF = —jp^(g 2 v + 3g 2 A ) dSn(S) -— [-f(-k)} x(-8 + k) F_ . (5.20) 

As discussed further in Appendix B, the T = part of this rate is obtained by taking 

-/(-*) -» i, 

rSr 0) = —£k- 3 (si + a&) f MS) f f x (-£ + 1) f_ . (5.21) 

7T Z7T r v 7 Jm J£+A fc 



VI. DISCUSSION 

We have derived expressions for the rate shifts, to order e 2 , in terms of thermal expec- 
tation values involving only electron field operators. The perturbation expansion of these 
expectation values has been carried out using the standard methods of thermal field theory. 
The "T = 0" terms have been carefully defined and removed. The remainder has been 
expressed in the form of convergent two-dimensional integrals. 

We gave formulae for only the neutron rates T n ; as we stressed earlier, if we had calculated 
the whole rate, rather than the finite temperature piece, we could use detailed balance 
directly to get a proton rate T p from T n . But since the present work is limited to the 
finite temperature part, we use the transformation A — > —A in the integrands to get the 
corrections to the proton rate from those for the neutron rate. We have performed the 

integrals in Eq's. Q5.9|), ( |5.10| ), and ( |5.13| ) to compute the total temperature-dependent part 



of the leading electromagnetic correction to the neutron rate as shown in Eq. ( |5.14j ). Their 
counterparts for the proton process have been computed in the same way after making the 
substitution A — > — A. We exhibit the results of the integration in the form of the fractional 
changes of the neutron rate AT^ /T^ and proton rate AT^ /T^\ where r^°) is the free 
neutron rate given in Eq. ( |2.30|) , while the free proton rate T-^) is obtained by the same 
calculation after the making the replacement A — > — A. The results are shown as a function 
of temperature in Fig. 1. 

We cannot compare these results directly with those in the literature, since previous 
authors have omitted the processes n + e + + u—>-p + '-f and p + 7— > n + e + + v. The first 
of these processes, since it has an outgoing photon, generates a T = part that has been 
subtracted as a part of the general T = subtractions in ( 5.1 4j) . We define a "revised" T part 



of the rate as that which must be added to the previous literature's T = part to get the 
complete answer. The "revised" T neutron rate is found by adding the contribution ( |5.21| ) to 



the the T rate plotted in Fig. 1. There is no T = piece of the reaction, p + 7— > n + e + + v, 
since here the photon is in the initial state, and therefore the proton "revised" T rate is the 
same as that plotted in Fig. 1. 

We plot the fractional "revised" rate corrections in Fig. 2. For comparison we also show 
data reconstructed from the figures of , shown as dashed curves; the latter are derived by 
combining the numbers shown in Fig. 11 of H, called the "finite T corrections", with those 
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of Fig. 16 of that reference, called the "finite temperature electron mass corrections." We 
see that at a temperature of 0.3 MeV our results nearly coincide with those of ||, but that 
the two families of results diverge rapidly at higher temperatures. We find positive rate 
corrections that are larger than those of Ref. for both processes, and we find much less 
difference between the neutron and proton corrections. Both of these changes lead to faster 
equilibration. 

Since the electromagnetic corrections are quite small, the detailed balance relation ( |2.9| ), 
expressed in terms of fractional corrections, is simply 

(Arf ) + Arf =°))/r(°) = (Arf) + Arf =°))/rf . (e.i) 

We can, in principle, try to check our answers with this statement by adding to our results 
the T = results of ||, shown in Fig. 8 of this reference, to our finite T corrections. 
Unfortunately this gets us into a spot of guesswork, since we cannot read the Arf 



=o) 



Arf =0 ) differences off of this plot quite precisely enough to get a clear answer. It would 
appear that these differences, as shown in the plot, are a little too small to compensate for 
the Arf) — Arf differences that we calculate and that are plotted in our Fig. 1. The same 
detailed balance test using the finite T results of 0, which shows a larger difference in the n 
and p corrections, appears to fail by more. We conclude that the approximations that were 
made in the calculation of the T = terms in Ref 's. Q and may deserve some further 
scrutiny. 

It is harder to make a detailed comparison of our results with those of Ref. [[?[] because 
of the way the latter are presented. However in Fig. 12 of Ref. |7j] we see large violations of 
the detailed balance requirement (|6.1|). 

We believe that with our formalism we have avoided the questions that have led to 
confusion and ambiguity in the literature. We should comment a little further on the wave- 
function renormalization questions that have plagued previous calculations. In the T matrix 
approach, with exclusive channels, the need for this renormalization arises from processes 
on an external line. Conceptually at least, one needs to sum the bubbles on the external line 
and calculate therefrom the residue of the energy shifted pole. The renormalization constant 
so obtained can then be expanded to second order in e. The complication in the case of finite 
temperature, with its preferred coordinate system, is that the multiplicative factor contains 
a transformation on the Dirac spinor indices, as worked out in [|T7]]. It is instructive to 
consider how these considerations would be embodied if we separated our result ( |4. 1 3|) into 



exclusive channels. The "external line" problem would present itself if we tried to evaluate 
the imaginary part of the square of the electron propagator, encountered in the third line 
of Eq. 13|) (with the omission of the spinor numerators) as Imfi^p) 2 — (E + ie) 2 ] -2 , by 
taking the imaginary part of one factor (a delta function) times the real part of the other 
(a denominator that is singular at the point of vanishing of the delta function). This is the 
structure that comes from taking a Feynman graph with a correction on an external line. 
We can see some of the benefits of the inclusive approach from the operations that follow 
our Eq. ( |4 . 1 3| ) , in which the representation of the imaginary part given by Eq. Q4.14| ) avoids 
entirely the ie delicacies required in treating exclusive channels. 

But there appears to be a deeper problem with the methods that have been used in the 
literature. In our Eq. ( ^4.13|) the propagators in question are the ordinary real time vacuum 
functions; the thermal effects have all been incorporated through the preceding formalism. 
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However, in j|] and || as well as in the subsequent papers on this subject, the individual 
graphs were calculated with real time thermal propagators of the form (for electrons) 



S(p) = (7P + my 1 + 2nn{\E\)(m - -?p)5(p 2 + m 2 ) . (6.2) 

When this propagator gets squared, it is truly meaningless; the interpretation is no longer 
just a question of following the ze's carefully. It appears to us that the square of the delta 
function has been implicitly treated in an arbitrary way in the treatments based on real 
time thermal propagators. This defect in the the real-time approach has been noted before 



in other contexts, for example in [TJ| and p0| . 

Finally, as we remarked earlier, the wave function renormalization is gauge dependent. 
In Feynman gauge, in which photons are emitted and absorbed by the proton, there is 
a temperature dependent part of the proton wave-function renormalization that must be 
calculated in order to get the correct answer in the exclusive approach. This appears to be 
missing in all of the published calculations. We conclude, from the above discussion, that 
it would be futile to try to make term by term comparisons of our results with the previous 
ones. However, we can say that most of the discrepancy comes from sources other than the 
new channel that we noted, 7 + p<->z/ + e + + n . This fact is made clear by Fig. 3, which 
shows the complete T = + T ^ contributions from these reactions. Since these terms 
obey detailed balance by themselves, there is no difference between the neutron and the 
proton fractional corrections shown in Fig. 3. 

A further use for the formalism developed in this paper would be in the calculation of the 
electromagnetic corrections to the charged current neutrino opacity in the environment of the 
supernova core. The neutrino reactions to be considered are the same ones that we have dealt 
with in the present paper. But the electrons are very degenerate, with the electron chemical 
potential in the 10's of MeV. Hardy [^T| has found, using some of the technology from 



the early universe calculations that we have discussed above, some evidence of significant 
electromagnetic effects in the rates. To get the complete answer to order e 2 one can adapt 
our equations to get the neutrino absorption rate in the presence of a non-vanishing electron 
chemical potential. 
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FIGURES 
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FIG. 1. The lower curve is (10 3 times) the ratio of the finite-temperature part of the leading 
electromagnetic corrections to the neutron disappearance rate Ar„ , given in Eq. ( 5.14] ) , to the 
free neutron rate Tn \ given in Eq. ( 2.30| ). The upper curve is the same ratio for the proton. 
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FIG. 2. Our results shown as solid lines compared to the results of Lopez and Turner [6] 
shown as dashed lines. As in Fig. 1, the curves describe the temperature-dependent part of the 
leading electromagnetic corrections to the rates divided by the corresponding uncorrected rates, 
Ar( T )/r(°). The upper curves refer to the proton rates, the lower curves to the neutron rates. 
As described in the text, we plot our "revised" result for the neutron rate obtained by adding 
the contribution (5.21) so as to account for the previously omitted process. To obtain the Lopez 
and Turner comparison figures, we combined their results for the "finite temperature radiative 
corrections", as shown in their Fig. 11, with their results for the "finite temperature electron- mass 
correction", as shown in their Fig. 16. 
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FIG. 3. Ratio of the previously omitted processes v + e + + n «-> j + p contribution to the 
neutron disappearance rate Ar n p, as given in Eq. ( |5.20 ), to the uncorrected rate T n ^ given in 
Eq. (|2.30| ), multiplied by 10 4 . Since the complete contributions of a new set of processes obey 
detailed balance, the corresponding ratio for a proton is the same. Note that the contribution of 
these processes is about an order-of-magnitude smaller than the Ar( r ) / Ar(°) ratios displayed in 
Fig. 2, and thus they cannot account for the discrepancies appearing in Fig. 2. 
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APPENDIX A: RATE FORMULA DERIVATION 



As discussed in the text, plasma modifications of the rates involve low-energy, soft pro- 
cesses. Thus they may be calculated using the limit in which the nucleons are treated as 
being very heavy and described by operators \Ed(r, £) and ^(r, t) that create and destroy a 
particle at the spatial point r at time t. The weak interaction Hamiltonian may be expressed 
as 



H w {t) =K(t)+K)(t) 



(Al) 



where 



G 



^ j (d 3 r){g v ijt(x) (1 - 75 )^e(x) ¥ n (x)* p {x) 

+g A ^(x)a, (1 - 75) Mx) ^{(x)a^ p (x)} 



(A2) 



Here ip e and ip v are the Dirac fields for the electron and (electron) neutrino, and ati = 7°7; 
is the usual Dirac matrix. 
The operator 



N u (t)= (d 6 r)^l(x)Mx 



(A3) 



measures the leptonic charge of the neutrinos, the number of v minus the number of v. The 
generic equal-time anticommutator for Fermi fields reads 



(A4) 



where the 5 function on the right implicitly includes the Dirac field indices. Using this 
anticommutator for the neutrino field gives the time derivative 



N„{t) = i[H w (t),N v (t)] 
K(t)-1C\t) 



(A5) 



The operators which appear here create and destroy electrons and neutrinos, and hence 
their expectation vanishes in the unperturbed plasma ensemble, which is diagonal in these 
particle numbers. The reaction rate T u appears in the additional linear response^ to the 
effect of the perturbation H w (t), 





11 dt 

' — 00 

dt 

-00 



N v (0),H w (t) 



/C(0)-/C t (0),/C(t) + /C t (t) 
Since K does not conserve particle number, 



(A6) 



7 We use here the general method advocated by Brown and Sawyer [22]. 
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<[/C(0) ,£(*)]> 



T 



T 



(A7) 



In virtue of the time transitional invariance of the ensemble, 



/C(0),/C f (t) 



/C(-t),/C f (0) 



(A8) 



Thus, we change the integration variable, t —>■ —t, use the antisymmetry of the commutator, 
and combine terms combine to obtain 



dt 



(A9) 



Although the expectation value that appears here describes a thermal ensemble of elec- 
trons and neutrinos, both of which are taken to have zero chemical potentials, we take 
the expectation value to have only a single nucleon, which we treat as a heavy particle lo- 
cated at the coordinate origin, r = 0. Thus the nucleon destruction operators in the weak 
Hamiltonian give a vanishing contribution unless their coordinate r = 0, and so the spatial 
coordinates of the electron and neutrino fields can be taken to vanish in the local interaction 
( |A2|) . Rotational invariance of the thermal ensemble implies that the (time-component) 
vector - (space-component) axial-vector interference terms vanish and that there can be no 
dependence upon the nucleon spin orientation. Hence the result is unchanged by averaging 
over the initial nucleon spin as well as performing the required sum over the final nucleon 
spin. Thus, we may replace the (direct product) 07 a m in the axial- vector - axial- vector 
contribution by 5i m . Moreover, as one can verify in detail from the structure of the neu- 
trino field thermal expectation values that will soon be presented, one can also replace the 
resulting (direct product) on on by 3. With these replacements made, both the vector and 
the axial- vector parts of the weak Hamiltonian involve the isospin lowering operator 

T_(t) = j ' (d 3 v)^ n (x)%(x) , (A10) 

and its Hermitian adjoint, the isospin raising operator 

T + (t) = Tl(t) = J (d 3 r)^ p (x)^ n (x) , (All) 

with 

^u = --G 2 w (g 2 v + ?>g 2 A ) J^dt 

( [^(0, t) (1 - 75) Ve(0, t)T_(t) , 4(0) (1 - 75) lM0)T + (0)] ) T ■ (A12) 

We define the real-time, generic free Fermi field thermal expectation values with vanishing 
chemical potential as 



(Mx)4(x'))P =\S {+) (x-x') 1 ° 



a/3 



(A13) 



and 



a/3 



(A14) 
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These function satisfy the free Dirac equation. Since these thermal expectation values involve 
Tre~P H • • •, with 



-I3H 



ip(r,t) =ip(r,t + if3)e 



(A15) 



the cyclic symmetry of the trace provides the boundary condition 



5W(r-r',t-07 Q 



a/3 



S { -\r-r',t + if3 -t')7° 



a/3 



(A16) 



This boundary condition, plus that provided by the equal-time anticommutation relation 
), determines the solution of the free Dirac equation to be 



S (±) (x-x')7° 



M_L P .P'M) 

(2vr) 3 2E(p) 
~E(j>) + 7 °m + a ■ pi n (=F#(p)) e- if ^-'') 



+ 



-E(p) — 7°m — a • p n (±E(p)) e 



-iE(p)(t-t>) 



(A17) 



Here 



and 



E(p) = Jp 2 + m? , 



n{E) 



1 



e^ + 1 



is the Fermi thermal distribution appropriate for an initial state, while 



n{-E) 



e PE + 1 



n{E) 



(A18) 



(A19) 



(A20) 



gives the Pauli blocking factor appropriate for final states. 

Since the rate ( |A12| ) already has the weak interaction taken into account, we may use 
the free field correlations ( [AIT] ) with m = for the neutrino fields. Since the electron 
thermal ensemble is parity invariant, a single factor of 75 does not contribute while t| = 1. 
Integrating over the neutrino solid angles and writing p = E u for the massless neutrino, we 
get 



T u = -Gl(gl, + 3gl) — j ^ dE v E 2 v n(E v ) j ^ dt 



,A > (2tt) 2 
-^Ve(0,t)T_(t)^(0)T + (0) 



Vl(0)T + (0)^ e (O,t)T_(t) 



(A21) 



The first term in the thermal expectation value only contributes when the initial nucleon 
state is a neutron which is changed into a proton by the action of T + (0); this term corresponds 
to the process u + n — > p + e~ plus the other processes related by crossing the leptons and/or 
including photon emission. These processes decrease the number v minus v in accord with 
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the overall minus sign out in front. Thus the positive rate T n for the disappearance of an 
initial neutron is given by 

dE v Eln{E v ) / dt (e + ^Ve(0,t)T_(t)^(0)T + (0) ) . (A22) 

-oo J — oo x ' 1 

Similarly, the rate for the disappearance of an initial proton is given by 

dE v E 2 u n{E v ) / dt (e- <B "Ve t (0)7+(0)^(0, t)T_{tj) (A23) 

-oo J — OO v ' 1 



APPENDIX B: ELECTRON SELF-ENERGY CONTRIBUTIONS 



The electron self-energy function for the thermal electron Green's function, is given, to 
one- loop order, by 



E(p) 



n=+oo 

-t £ 

n=~ oo 



(d 3 k) m - j{p + k) 2 
(2 7 r)3 7 ' 1 (p + A;) 2 + m 2 ^ 6 M j ' 



(Bl) 



The four vectors here are in Euclidean space with the scalar product 7(2? + k) = 7 • (p + 
k) + 74 (p + fc) 4 , 74 = — 27 , and k A = uj n = 2nnT while the external energy takes on the 
values p 4 = (2m + l)nT. The radiation gauge is chosen where 



Di m (k) 



Sim. — 



kik m 



fc2 ' 



(B2) 



and 



D 44 (k) = p . (B3) 
The frequency sum is conveniently performed by expressing it as a contour integral, 

where the contour C runs from —00 to +00 just below the real axis and then returns to —00 
just above the real axis. For the part of the contour just below the real axis, we write 



1 



1 



oi/3uj _ Y ' 



(B5) 



with the second term exponentially damped in the lower-half plane where we shall close the 
contour for it. For the upper-half portion of the contour, we write 
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and close the contour for the second term in the upper-half plane where it is exponentially 
damped. In this way, we achieve the separation 



E(p) = ^°\p) + ^ T \p), 



(B7) 



where 



£(°)(p) 



dw m - 7 (p+fc) 

2nJ (27T)^(p + ky + m^ u ,lu[ ) 



(B8) 



is the vacuum self-energy function (in Euclidean space). This vacuum part must be regulated 
and renormalized. The remaining, temperature-dependent piece £( T )(j») is a finite, well- 
defined function which vanishes when the temperature vanishes. 

In closing the contours for the thermal part of the self energy, poles are encountered at 
uj = ±ik, where k = |k|, producing a factor of the Bose distribution 



/(*) 



1 



ol3k _ I ■ 



(B9) 



Poles are also encountered at uj — ±iE(p + k) — p 4 , where E(p + k) = \J (p + k) 2 + m 2 . 
Since (3p± = (2m + l)n, with exp{iPp 4 } = —1, these poles are accompanied by a factor of 
the Fermi distribution 



n(£(p + k)) 



e /3£(p+k) + l ■ 



(BIO) 



Separating out the Coulomb and radiative pieces, and continuing to real energy by setting 
P4 = —iE, we obtain 



where 



4 T) (p,£) 



7/ 



(d 3 k) 
~{2nf 



2k 



f(k) 



m 



£^p) = 4 t V) + 4 t) (p), 



(film klkrnj 

7 • (p + k) + 7°(E + k) m - 7 • (p + k) + 7 °(£ - k) 



(Bll) 



+ 



E(p + k) 2 - (E + k) 2 E(p + k) 2 - (E - k) 2 

~m - 7 • (p + k) + 7°£(p + k) 



2E(p + k) 



n(£(p + k)) 



k 2 - (E — E(p + k)) 2 



and 



p,E) = -2e 2 



m - 7 • (p + k) - 7°£(p + k) 
k 2 — (E + E(p + k)) 2 

(d 3 k) 



1 



(2vr) 3 2£(p + k) 



n(£(p + k)) 



m + 7 • (p + k) 

k 2 



(B12) 



(B13) 
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A partial fraction decomposition in the radiative piece produces a convenient alternative 
form: 



i T) « 



ll 



+ 



(rf 3 k) 

V lTn '" r " m ) 4£;£(p + k) 
m - 7 ■ (P + k) + l°E(p + k) 

7(fc)+n(£(p + k)) /(fc)-w(£(p + k)) 
_ £(p + k) - E - k E(p + k)-E + k 

m - 7 • (p + k) - 7°£(P + k) 

7(A;)-n(£(p + k)) , f(k) + n(E(p + k)) 



(B14) 



£(p + k) + £ + fc £(p + k) + E - k 

Thus we find that the gamma-matrix trace in the basic rate formula ( 4.15|) involves, for the 
radiative part, 



4 (Sim - kkm) tr (m - 7p)7°(m - 7p) 

m - 7 • (P + k) ± 7°-E(p + k) 



7/ 



7r, 



= —AEm 2 + Ek 2 — (E/k 2 ) [E 2 (p + k) - £ 2 (p) 
±2£(p + k) [E 2 + E 2 (p)" . 
The trace for the Coulomb part is given by 

tr (m — 7p)7°(m — 7p) [m + 7 • (p + k)] 
£ 2 (p) + p-k 



(B15) 



1 

4 
2£ 



(B16) 



We begin by calculating the ImS contribution in the rate formula ([4.151) , which we 
shall call r^_ e . The complete self-energy function (with no T = subtraction) enters here. 
Returning for a moment to recall the discussion of our evaluation of the frequency sum (|B4Q . 
it is easy to see that the complete S function can be obtained from the result (|B14j) for E^ T ) 
by making the replacements, f(k) — ► [f(k) + 1/2] and n(E(p + k)) — > [n(E(p + k)) — 1/2]. 
We incorporate this replacement and insert the imaginary part of the result ( [B14| ) evaluated 
with the trace formula (B15| ) into the basic rate formula (|4.15| ). We interchange the p and 
k integrals and then replace the p integration variable by q = p + k. In this way, we obtain 



^(7) 



p 2 n 2 



, . r+oo 

[9v + Z9a) dE X (E) 
x ' J— 00 



(27T) 2 

(d 3 k) 1 r (d 3 q) 

(2tt) 3 k J (2vr) 3 E(q) [E 2 (q - k) - E 2 ] 2 



2tt 



[f(k) + n{E{q))\ [r + S{E + k- E{q)) -r_5{E-k + E{q))} 



+ [f(k) - n{E{q)) + 1] [r+ 5 {E - k - E{q)) -r_5(E + k + E{q))} 



(B17) 
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The only angular dependence that appears here is in 



E{q-k) 2 = E{q) 2 + k 2 + q-k, 



(B18) 



which enters in the denominator and in the definition ( P15|) of r±. Taking k as the z axis 
and performing the solid angle integral for q and then the remaining solid angle integral 
for k, and using the delta functions to remove the E integral gives, after one final variable 
change to £ = \J q 1 + m 2 , and some algebraic effort, 



r (7) o e 

- 1 - ne-e " 



(g 2 v + %g 



, oo rlh roo 



(2tt) 5 V y ' ~ aA ' Jo k 
[f(k) + n{£)\ [ X (£ — k) + x(-£ + k)] F_ 



+ [f(k) - n{£) + 1] [ X (£ + k) + x (-£ - k)} E 



+ I > 



(B19) 



where 



with 



and 



F± = A±kB 



A 



2£ 2 + k 2 



B = 2£ln (^)- 4 «- 



(B20) 



(B21) 



(B22) 



We have also made use the function defined in Eq. ( |2.29| ), which we repeat here for conve- 
nience, 



X{E) = n(E - A)n{-E){E - A) 2 . 



(B23) 



To facilitate the removal of the "T = 0" piece and to also put the result in a form that 
clarifies its structure, we note that 



and 



f(k) + n{£) = f(k) n{-£) n{-£ + k)~ l 
= -f{-k)n{£)n{£ -k)' 1 



f(k) - n{£) + 1 = -f(-k) n{-£) n{-£ - k)~ l 
= f(k)n(£)n(£ + k)- 1 . 



(B24) 



(B25) 



Here 



/(-*) = ! + /(*) 



(B26) 
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and 



n(-£) = l-n(£) (B27) 

describe the occupation factors appropriate for final states: l + f(k) gives the Bose enhance- 
ment factor and 1 — n(£) the Pauli blocking factor. Using these identities, we obtain 

| [f(k) n{-£) x(S - k) - f(-k) n{£) x(-£ + k)] F_ 
+ [-f(-k) n(-S) x(£ + k) + f(k) n{£) x(-£ - k)] F+\ , (B28) 



where 



X(E) = n(E - A) (E - A) 2 . (B29) 



The energy integration variable E in Eq. ( |B17| ) differs from the neutrino energy by a positive 
constant, E = E v + A. Neutrinos in the initial state are described by positive values of the 
neutrino energy, E v > 0. Thus the terms in Eq. ( |B17| ) involving 5(E — k +•••), where 
k — |k| is an on mass shell, real photon energy, correspond to the production of photons in 
the final state. These real photon contributions appear in Eq. ( |B28| ) in the terms involving 
x(— £ + k) and x(£ + k). These are just the terms that have the occupancy factors f(—k) 
that are appropriate for produced photons. It is a simple matter to check that the Fermi 
factors n(±E) are just those needed in the various photon absorption and emission processes 
described by Eq. (|B28|) , but we shall not bother to enumerate then here. We just note that 



the "T = 0" part of the rate is obtained by setting f(k) = 0. Subtracting this part is 
equivalent to the substitution —f(—k) — ► f{k), and so the T ^ 0, photon thermal bath 
contribution is given by 



p(7.7) o 
1 ne—e ^ 



[n(-£)x(£ ~k)+ n(£)x(-£ + k)] F_ + [n{-£)x{£ + k) + n{£)x{-£ - k)] E 

m 



(B30) 

Although this set of terms has no ultraviolet divergences, since f(k) k~ l — > k~ 2 as k — > 0, it is 
infrared divergent. This contribution is rendered finite in the infrared by other contributions 
that we are about to compute. 

We pause for a moment from our development to identify the contribution from the rate 
e + + v + n — >7 + p which has been omitted from all previous papers. Since it involves real 
photon production, it must involve the terms containing —f(—k) in Eq. ( |B28| ). The second 



term involving the + function describes an incident neutrino of energy E v = £ + k — A 
corresponding to a produced electron of energy £ rather than an initial positron of this 
energy. Hence it is the first term containing —f(—k) that contributes to the process in 
question, with the photon energy starting at k — £ + A corresponding to an incident 
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neutrino with zero energy. Thus the previously forgotten contribution to the rate is given 
by 

r BP = 2^((£ + 3(&) J m dSn(S) J £+aT [-f(-k)}x(S + k)F.. (B31) 
The T = part of this rate is obtained by taking —f(—k) — > 1, 



<T=0) 
nF 



(2tt) e 



(<7y + 3p 



d£n(£) 



x(-f + fc) F_ 



(B32) 



We turn now to evaluate the on mass shell E = ±E(p) contributions that appear in 
the last two lines of the rate formula (|4.15 ). This piece of the rate formula involves the 
temperature-dependent self-energy function in the form 



E (T) (p, E) = tr (m - 7p) 7 (m - 7p) S (T) (p, £) . 



(B33) 



We separate out the parts that involve the photon and electron phase-space densities by 
writing 

m (p, E) — Ef ) (p, E) + Ef ) (p, £) . (B34) 
We use of the previous results QB13Q - QB16Q and a little algebra to write the two parts as 



(2vr) 3 jfc 



m 2 + p • k + 



+ 



(p ■ kg 

k 2 
1 



£ 2 + £ 2 (p) 



£ 2 (q) - (E + k) 2 E 2 (q)-(E-k) 2 
E+k E-k 



+ 



E 2 (q) - (E + k) 2 E 2 (q)-(E-k) 



(B35) 



and 



Ef)(p,£) = -2e 2 



(rf 3 q) 1 



(2tt)3 £(q) 



f 4£ 




\ fc 2 " 


£ 2 (p) + p-k 



A; 2 - 4m 2 - — 
1 



£ 2 (q) - £ 2 (p) 



+ 



+ 2£(q) £ 2 + £ 2 (p) 



(£(q) - E) 2 - k 2 (E(q) + E) 2 - k 2 
1 



(£(q) - £) 2 - fc 2 (£(q) + £) 2 - A; 2 



(B36) 



Here we have again introduced the variable q = p + k and used it as the integration 
variable for the electron contribution to the self energy ( |B36|) . We have used the different 
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integration variables k and q for the photon and electron contributions because then the 
angular integrations appear in a simple form. The terms on the first line on the right-hand 
side of Eq. ( |B36| ) come from the Coulomb contribution ( |B13| ) to the self energy. 

To proceed with our computation, we note that the self-energy function evaluated on the 
mass shell is gauge invariant^ and gives an energy shift .0 For positive energies (particles) 
with spinors that obey 



V m a (jo) u\(p) = m->yp, 



we have 



(B37) 



(B38) 



2EAE = u x (p)Y>u x (p), 
while for negative energies (antiparticles) with spinors 

Y, x v \{p) v\(p) = m + 7P , (B39) 

the energy shift is 

2EAE = v x (p)J:v x (p), (B40) 
Thus, in either case, one of the pieces of the rate formula Q4.15Q may be expressed as 



8 E 2 AE {T) (p) = tr 7° (m - jp) £ (T) (p, E) (m - jp) 
= ^ T \p,E) 



E=±E(p) 



(B41) 



E=±E(p) 



Accordingly, we write the last, on-mass-shell contribution to the rate formula ( [4.1 5| ) as 

rlft-l + rie_e, where 



ri A e 2 = Gl (gl + 3<?i) \ £j dE J 



(d 3 p) 
(2tt) 3 







[S (E - J5(p)) + 5 (E + E(p))} AE^ T \p) — X (E) 



(B42) 



and 



r>2 
(2vr) 2 



(g 2 v + 3<? 



— oo 



(2tt) 3 2E(p) 
d 



[5 (E - J5(p)) -6{E + J5(p))] X (E) ^ [J5T 1 £^(p, E) 



(B43) 



8 That is, the result would not be altered if we changed our use of the radiation (or Coulomb) 
gauge to a relativistic gauge such as Landau of Feynman gauge. 



See, for example, the discussion of Eq. (8.7.13) in Section 8 of |23]. 
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Referring back to the form of the uncorrected neutron rate ( |2.30|) , we see that Eq. ( |B42|) 



describes precisely the correction to that rate which results from the (gauge-invariant) energy 
shift AE( T \p). As we shall see, this temperature-dependent correction is well defined 
with no infrared divergences. The remaining contribution (|B43|) , which has the form of a 
temperature-dependent wave function renormalization, is not gauge invariant, and it also 
suffers from an infrared divergence. This contribution by itself has no physical meaning. 
It combines with previous contributions to form a gauge-invariant result with no infrared 
divergence. 

In parallel with the partition ( |B34j ) , we write 

AE^ (p) = AM T ) (p) + AEjP (p) , (B44) 



with the total energy shift given by Eq. (B41|) . Using Eq. ( |B35|) for the photon part gives, 
after some calculation, 

^) = i/wl /(l) ' (B45) 

where the positive/negative energy has the mass shell values E = ±E(p). Similarly, some 
work with the electron part ( |B36|) yields 



2 9 

m m 



[E(q) + E(p)} 2 -k* + [£(q) - E(p)} 2 - A? J ' ^ ' ° ' 



Averaging over the solid angle of q produces 

These energy shift corrections are perfectly finite: They have neither ultraviolet nor infrared 
divergences. They also vanish in the zero temperature limit. Moreover, as remarked before, 
the energy shifts are gauge invariant. This gauge invariance may be confirmed by explicitly 

frp\ 

computing the self-energy function in the Feynman gauge, Y> F (p, E). This function is ob- 
tained by making the replacement 5i m —kik m — > in Eq. ( |B12| ), where <? M „ is the Minkowski 
four- dimensional space-time metric with the spatial gamma matrices ji, 7 m replaced by the 
space-time matrices 7 M , j u . The Coulomb contribution ( |B13| ) is now, of course, omitted. It 
is a simple matter to evaluate this replacement. Passing to the mass shell, —p^p^ = m 2 , 
with effectively 7 M — > p^/m, it is easy to show that 

Z { P(p,E)^-AE( T \p), (B48) 
m 



where AE^(p) is given by Eq's. (P44 ) - (P46 ). (Our previous work shows that E/m is 



the correct factor to relate the energy shift to the mass-shell self energy.) Gauge invariance 
is the reason behind the cancelation of the 1/k 2 terms when the electron contribution ( B36|) 
to the self energy is evaluated on the mass shell. The resulting energy shifts are the shifts in 
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the positions of the poles of the thermal electron Green's functions — they give the change 
in the energy needed to create an electron or positron of momentum p in the plasma relative 
to the energy needed to create them in the vacuum. 

The energy shifts are physical quantities. Thus we let them stand as separate corrections 
and turn to the temperature-dependent "wave function renormalization" contribution (|B43|) 



which is gauge-dependent and which is not infrared finite, a contribution that must be 
combined with the previous result to obtain a physically relevant contribution to the neutron 
rate. Some calculation with the photon part (|B35|) of the self energy yields 



_d_ 

dE 



2e 2 1 [™ dk 

t 2 / -rf( k ) A - 

7T p JO K 



(B49) 



E=±E(p) 

where A is the function defined before in Eq. QB21J ) but with different labels, 

/ E(p) +p' 



A 



2E(py + k z 



In 



\E{p) - p ; 



-ApE(p), 



(B50) 



The electron part (|B36|) of the self energy yields 



_d_ 

dE 



E-^P( P ,E) 



2e 



E=±E(p) 

, E\q) 
pq 



oo 



tt 2 ./., dq E(q) 



n{E{q)) 



E{p) 



E 2 (q)-E 2 (p) 



In 



p + q 



E(q) E\q)+EZ{p) l J[E{p)E{q)+_ 



l 2 — m 4 ' 



E(p) 



2pq 



\ [E(p)E(q) — pq] — m 4 / 



(B51) 



For the photon contribution ( |B49|) to the wave function term (|B43| ) , we integrate over 
E and make the simple notational change p — > q with, as before, £ = E(q). Then, since 
qdq = SdS and, integrating over solid angle, 



(rf 3 P ) 71 

(2tt) 3 2E(p) 



1 

47T 



q dS 



(B52) 



Thus, the photon contribution to the wave function correction (|B43|) is 
TZte = 2^(g 2 v + 3g 2 A ) jT ™f(k) £ d£ A {-2 X (£) - 2 x (-£)} 

= 2 T2#(^ + 3 ^) I T f{k) L d £ A i-M-£)m-M£)x(-£)} ■ (B53) 

As discussed in the text, this adds simply with the real photon contribution Eq. ( |B30| ) to 
produce an infrared finite result. The second wave function contribution comes from placing 
Eq. ( |B5ip in Eq. ( |B43j ) • Removing the first energy integral by the delta functions, performing 
the angular integral, and changing integration variables from momentum to energy with the 
notational change q — > p', E(q) — > E', gives the result 
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£ 2 G\y 



(27T) 



(g 2 v + 3<? 



£' 2 - _E 2 
~~E 



App' 



dE 



2 £' 2 In 



dE'n(E') [ x (E)+x(-E) 
p + p \ 



E' 2 + E 2 



In 



p — p' J 
[EE' + pp'} 2 — m 4 
[EE' — pp'] 2 — m 4 



(B54) 



This piece by itself is not gauge invariant. As discussed in the text, adding it to Coulomb 
interaction between the electron and proton, the rate contribution (|3.9|), produces a well- 
defined, gauge-invariant contribution. 
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